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We investigate the electromechanical coupling in 2d materials. For non-Bravais lattices, we find 
important corrections to the standard macroscopic strain - microscopic atomic-displacement theory. 
We put forward a general and systematic approach to calculate strain-displacement relations for 
several classes of 2d materials. We apply our findings to graphene as a study case, by combining 
a tight binding and a valence force-field model to calculate electronic and mechanical properties 
of graphene nanoribbons under strain. The results show good agreement with the predictions of 
the Dirac equation coupled to continuum mechanics. For this long wave-limit effective theory, we 
find that the strain-displacement relations lead to a renormalization correction to the strain-induced 
pseudo-magnetic fields. Implications for nanomechanical properties and electromechanical coupling 
in 2d materials are discussed. 


I. INTRODUCTION 

Linear continuum elasticity provides a valuable ba¬ 
sis for the investigation of the mechanical properties of 
atomic monolayer materials [U [2]. Elastic theory is also a 
key element for understanding the relation between the 
material deformations and the corresponding modifica¬ 
tions of its electronic structure[^-[l]. Roughly speaking, 
strain changes the interatomic distances, thereby mod¬ 
ulating the overlap of electronic orbitals of neighboring 
atoms and modifying the electronic properties of the ma¬ 
terial. To model and engineer this electromechanical cou¬ 
pling it is necessary to correctly relate experimentally 
controllable macroscopic deformations, parametrized by 
the strain tensor u, to microscopic atomic displacements. 

We find that for materials with a crystal structure with 
a basis there is an important correction to the standard 
strain-displacement relations. We show that by applying 
strain to a 2d monolayer material, the nearest-neighbor 
vectors connecting the atoms transform as r^- —>• (fla -|- 
u) • Vij -\- A, where is the n x n identity matrix and 
A is a vector that depends on the material deformation 
energy. 

In this Letter we provide a transparent and systematic 
approach to calculate the strain-deformation relations for 
any kind of 2d material, provided that the deformation 
energy can be parametrized in terms of the bond vectors. 
We show the significance of our findings by contrasting 
the elastic bulk properties of graphene and black phos¬ 
phorous calculated with and without the proposed cor¬ 
rections. 

We study the effect of the strain-displacement correc¬ 
tions on electronic degrees of freedom in two applica¬ 
tions. In the first one we incorporate the modified strain- 
displacement relations in the extensively used k • p the¬ 
ory for graphene, where strain can be represented by an 
effective gauge field [^]7]. We show that a correct as¬ 
sessment of A renormalizes the strength of the gauge 
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vector-potential by a factor 3/5, which means a factor 
(3/5)^ for scattering rates relevant to transport prop¬ 
erties 0 0. A similar correction has previously been 
identified in the study of the elastic properties of car¬ 
bon nanotubes, where the strain-deformation relations 
were inferred from the analysis of the acoustic phonon 
modes [51 [10]. In the second application we study the ef¬ 
fect of our findings on the gap engineering of graphene 
nanoribbons using strain. 

In summary, our study shows that quantitative esti¬ 
mates of the electromechanical coupling in 2d materials 
- as is often required in strain-engineering applications - 
requires careful consideration of the lattice deformation 
on a microscopic level. 

II. RESULTS 

The atomic lattices of 2d materials are characterized 
by primitive unit cells (PUCs) with a set of basis atoms 
and two primitive lattice vectors ai and a 2 . In a homo¬ 
geneously strained sample it is sufficient to consider the 
deformation of a single unit cell under an applied strain. 
The deformed primitive lattice vectors read 

a' = (Us -I- u) • ai, (1) 

where u is the strain tensor and i = 1,2. This equation 
also holds for nonuniform strain, for which u = u(r) varies 
on length scales much longer than the lattice parameter. 
In general, the lattice vectors can be expressed by a linear 
combination of the bond vectors Vjk, 

^ ^ CijkVji^ , ( 2 ) 

rjfce{rD 

where {rd is the set of bond vectors contained within 
each unit cell and Cijk is a tensor encoding the lattice 
connectivity. We address the case where the strain field 
deforms the lattice, but preserves its connectivity. We 
show that for materials with a non-Bravais lattice struc¬ 
ture, the modifications of the bond vectors due to strain 
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do not follow Eq. Q . Let us start from the most general 
relation 

Ty- = (^3 + U) • Tij + Ay , (3) 

where the vectors Ay contain the differences between the 
strain displacement relations for and ry. The number 
of independent vectors Ay- is determined by the number 
of basis atoms and ^efr} = 0. For a Bravais 

lattice Ay- = 0. Each additional atom introduces three 
degrees of freedom, that can be expressed by an in-plane 
vector Aj| and an out-of-plane component Aj^ z. The 
number of independent displacement vectors obtained in 
this way can, in general, be reduced by enforcing the 
lattice symmetry and connectivity. To be specific, we 
consider the most important 2 d materials [llj under cur¬ 
rent investigation, divided in three groups, shown in Fig. 

EH 

The first class contains graphene and graphene-like 
materials, such as boron nitride, as well as non-planar 
materials, such as silicene and germanene. The PUC 
consists of two basis atoms and two lattice vectors ai = 
ao (3/2, •\/3/2,0) and a 2 = oq (3/2, —•\/3/2, O), see Fig. 
©a). The three nearest neighbor vectors are given by 
ri = (—oo, 0 , Q) — hz, r 2 = ai-|-ri and ra = a 2 -|-ri, where 
h is the PUC corrugation height. The lattice vectors are 
written in terms of the bond vectors as ai = r 2 — ri 
and a 2 = ra — ri. The bond vectors transform as 
r' = Ti -f u • Ti -f Ai, with the constraints 0 = A 2 — Ai 
and 0 = A 3 — Ai. Hence, A^ s A. The latter is conve¬ 
niently expressed by in-plane deformations Ay and the 
corrugation change A_l 2 . Graphene is a special case for 
which Ax = 0 . 

The second class contains the transition metal 
dichalcogenides (TMDCs). Here the PUC contains three 
atoms, one transition metal (Mo or W) and two chalco- 
gens (S or Se), see Fig. [S^b). The transition metal 
has six nearest neighbors. The bond vectors are la¬ 
beled r 2 ^u/d and ra ,^/^^ where the subscript u 

or d indicates whether the bond connects to the up¬ 
per iu) or lower (c?) layer of chalcogen atoms. The 
lattice vectors are given by ai = Y 2 ^u/d — '^i,u/d and 
a -2 = '^ 3 ,u/d ~ ^i.u/di- The bond vectors transform as 
^'i,u/d = U,«/d + ii'U + ±Axx/2z. The vectors A,j| 
now satisfy 0 = A 2 ,|| — Aij and 0 = Aaj — Aiy|. Conse¬ 
quently Aj II = A|j. By inversion symmetry, vectors con¬ 
necting to the upper layer and the lower layer transform 
in the same way, but with opposite signs in the out-of- 
plane direction. Hence, the strain-displacement response 
of TMDCs is also characterized by a single in-plane vec¬ 
tor A|| and the change in inter-chalcogen distance Ax. 

The third class consists of a single material, phos- 
phorene, an atomically thin puckered material derived 
from layered black phosphorus. Its PUC contains four 
atoms and six unique bond vectors, see Fig. [sT|(c). Four 
bond vectors are situated in the puckers and are denoted 
by Vi u/d with z = 1, 2 and subscripts u/d. The two 
remaining vectors connect the upper and lower puck¬ 
ers, and are labeled Vcc and Vcc'- We find that ai = 


-Tec + ri,„ -k Tec' + r 2 ,d, and a 2 = -r 2 ,„ -k ri,^. Al¬ 
ternatively, we also write ai = —Ycc + i’ 2 .« + r/c + ^i,d 
and a 2 = —r 2 ^d + ^i,d- Hence, we obtain the constraints 
0 = — Acc + Ai^u -k Acc' + A 2 ,d) 0 = —r 2 ,u -k ri_„. This 
leads to four equations for the vectors A^. Using in¬ 
version symmetry and the same arguments as for the 
TMDCs, the strain displacement relation in phosphorene 
is also characterized by a single in-plane vector Ay and 
a change in vertical inter-pucker distance Ax. 

We conclude that upon application of strain to all these 
2d materials, due to lattice connectivity and symmetry 
constrains, the strain-displacement correction is simply 

Ay- = A|| + Axz, (4) 

independent of i and j. Both Ay and Ax depend on the 
applied strain and on the interatomic interactions. Their 
calculation requires a microscopic model to account for 
the deformation energy. For uniform strain, one can de¬ 
termine the strain-displacement correction using the geo¬ 
metric considerations of the previous paragraph and, for 
instance, first principles calculations. For non-uniform 
strain and/or finite lattices, this procedure becomes com¬ 
putationally prohibitive, and one has to resort to semi- 
empirical force-field models. 

In the following, we put forward a systematic and 
transparent approach to calculate Ay and Ax for any 
2d material whose deformation energy is described by a 
force-field model parameterized by the bond vectors. Va¬ 
lence force models (VFMs) are a convenient choice since 
they offer good accuracy at low computational costs. 
There are various VFMs for the different 2d materials 
[IMl]. Let us focus on graphene-like materials, which 
can be accurately addressed using the VFM introduced 
by Perebeinos and Tersoff[T4] to describe the interactions 
between sp^-bonded carbon atoms. The deformation en¬ 
ergy is given by mi 

^def = 4 E 

dr) 

^ ijGi iJ<kGi 

+ — (Sr.j) (Sci ki) 

+ ^ V ( V 

where = |rC| is the bond length, fry = rC — oq 
is the change in bond length and 6cijk is defined as 
Scijk = cos9[ jk — cos 9ijk- Here 9ijk is the angle be¬ 
tween atoms z, j and k with atom z as apex in equilib¬ 
rium, while 9'- jf. is the angle in the deformed lattice. The 
summations in Eq. ([^ follow the convention: (a) j S z 
indicates that the index j runs over the three neighbors 
of atom z; (b) for j < k G i both j and k are neighbors 
of the atom z, and are ordered to avoid double counting. 
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FIG. 1. Lattice structures of the three classes of 2d materials considered, a) The unit cell of graphene-like materials are defined 
by two lattice vectors ai and sl 2 and three nearest neighbor vectors ri, r 2 and rs. The material can in general be corrugated, 
signified by h in the lower figure. Graphene is contained as the case h = 0. b) The unit cell of TMDGs consist of two lattice 
vectors and three basis atoms, which define six nearest neighbor vectors. Due to the layered structure of the material, it has 
an effective thickness h. c) The phosphorene lattice is characterized by two lattice vectors and four basis atoms, defining ten 
nearest neighbor vectors. Its puckered structure gives rise to the effective thickness h. 


leaving three possible terms; {c) j k < I G i leaves three 
terms for each i. The first and second terms in Eq. ([^ 
give the energy cost of stretching and bending bonds as 
in a Keating model[T7]. The third and the fourth terms 
couple stretching of different bonds and couple stretch¬ 
ing and bending, respectively. The last term is related to 
out-of-plane displacements. In Ref. m there is an addi¬ 
tional term which penalizes misalignments of neighboring 
TT-orbitals that we neglect here. 

For a given strain, we minimize Kdef with respect to 
Aj| and Ayz to obtain 


aoK{h) 

2 



Ax = ± aoK±{h){uxx + Uyy) , 


(6a) 

(6b) 


where k(/i) and K±{h) are functions of the unit cell cor¬ 
rugation height h and characterize the changes in the 
PUC internal structure, in addition to the trivial shift of 
atomic positions when the PUC is strained. To lowest 
order in h/ao, 


n{h) 

n±(h) 


9/3e - 4/3,. + 2/3.2 

9/3c + 4/3r ~ + ^Prc) 

h 2/3r + 2/3^2 + 3/3rc 
ao 9/3y 


(7a) 

(7b) 


Upon elimination of Ay and Ax, the microscopic de¬ 
formation energy per unit area, £ = Edei/A, becomes 
the standard expression for the continuum elastic energy- 
density of an isotropic membrane [3], 


£ = —X{uxx + Uyy) ^{Uxx T Uyy -j- 2Uxy) 

Y^2d 

= ^^^{Uxx T Uyy 2vUxxUyy -f (1 — v)Uxy), (8) 


reconciling our results with previous works on the long- 
wavelength elastic behavior of graphene [T1 [TOl . Here, 
A and /r are the Lame parameters of the material, and 
Y 2 d and v are the 2d Young modulus and Poisson ratio, 
respectively. They are related via 


Kgd — A -|- 2/i , 


and 


A 

A -t- 2fjj 


(9) 


We find that the Lame parameters, given in terms of the 
microscopic parameters of Eq. ([^, are 


A = 


•\/3q 


8^^ - ISMr + 4/3, {13x2 - 3/3.c) 
9/3c + 4/3j. — 2 (/3r2 + 3/3rc) 


9 / 32 , - mMr2 - 4 / 3^2 - l2l3rcPr2 


3^3 

% 


9/3c -|- 4/3^ — 2 (/3r2 + 8l3rc) 
- 2/3c/3r2 - lire 


9jic + 4/3r — 2 IPr2 + ^Pre) _ 


(10a) 

(10b) 


For graphene, the parameter values reported in Ref. [T3I 
together with Eqs. ,0 lead to K = 0.39 and kj_ = 0. 
Moreover, Eqs. ( [I0| give A = 4.4 eV/A^ and fj, = 8.8 
eV/A^, which implies Y 2 d = 352 N/m and jz = 0.2. 

Neglecting the bond-bond and bond-angle correlations 
in Eq. ([^ (setting Pr 2 = Pro = 0), we find that K,{h) 
depends only on the Poisson ratio, n{h) « 2jz/(1 -|- v). 
Consequently, k( 0) vanishes for materials with a negli¬ 
gible Poisson ratio (where /i A), and is limited from 
above by k( 0) < 2/3 for isotropic materials (for which 
[3] iz < 1/2). In this simplified VFM, k( 0) agrees with 
the findings of Ref. El and may serve as a rough esti¬ 
mate, since it can be readily obtained form the Pois¬ 
son ratio which is a macroscopically observable quantity. 
For V = 0.2 this approximation gives n « 1/3, which 
agrees reasonably well with the value obtained from the 
full set of parameters (k = 0.39). Similarly, we find 
Kx(h) « -2{h/ao){/3r/{9Pr)). 
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FIG. 2. Structure of an armchair nanoribbon. The high¬ 
lighted area represents the system unit cell. The ribbon con¬ 
sists of Ni X N 2 cells, where Ni counts the cells in ^-direction 
and N 2 in j/-direction. In the lateral direction periodic bound¬ 
ary conditions are imposed at the edges. 



applied strain Uxx 


The approximation k « -|- v) is expected to 

hold also for TMDCs, due to the hexagonal struc¬ 
ture of the lattice when projected onto the monolayer 
plane. For example, using ;/moS 2 = 0-27 [H], we obtain 
«^MoS 2 ~ 0.43, which is slightly larger than the estimate 
for graphene. For phosphorene, the strain-displacement 
relation no longer obeys the simple form of Eq. §, 
since A shows a directionality reflecting the material 
anisotropy. Nonetheless, we can still use our theory to 
estimate the elastic properties of phosphorene (see Sup¬ 
plementary Material). Using a VFM developed for lay¬ 
ered black phosphorus j2], we obtain sound velocities 
Vxx = 3508 m/s, Vyy = 8147 m/s, and Vxy = 3707 m/s in 
good agreement with Ref. [5] (without correction Vxx and 
Vxy are overestimated by a factor 2). The Young moduli 
are Fj, = 17 N/m and Yy = 94 N/m (without corrections 
Yx = 58 N/m and Yy = 95 N/m). These values are in 
good agreement with DFT calculations I3 0IS], where 
Yx = 24... 29 N/m and Yy = 88... 102 N/m. Hence, 
a correct treatment of the strain-displacement relations 
is necessary to account for the predicted anisotropy of 
phosphorene [5^ . 

We illustrate our results by studying armchair 
graphene nanoribbons (AGNRs) stretched along the lon¬ 
gitudinal direction. See Fig. for a sketch of the setup. 
The equilibrium conhguration is obtained by minimiz¬ 
ing Edef, Eq. (§. The mechanical energy density as 
a function of the applied strain Uxx is shown in Fig. 
[^a). Numerical results are denoted by symbols, whereas 
full and dashed lines correspond to analytical results ob¬ 
tained from Eq. (|^ for k = 0.39 and k = 0. By setting 
At = 0, one obtains a slightly larger value of Y 2 d such that 
the stretching energy is overestimated. We also calculate 
the resulting strain Uyy in the transversal direction. As 
the ribbon contracts upon stretching, Uyy is negative for 
Uxx > 0 and decreases with increasing strain. This is 
shown in Fig. [^b). The slope of the curve at Uxx = 0 
gives the Poisson ratio, which agrees very well with the 
estimate obtained from Eq. The standard approach, 
for which k = 0, gives za(k = 0) « 0.16, thus underesti¬ 
mating the Poisson ratio. 



applied strain Uxx 


FIG. 3. Strain Uxx dependence of (a) the elastic energy and 
(b) the resulting strain Uyy for AGNRs of different widths 
(Ai = 300 in all cases). Symbols denote numerical results. 
The full (dashed) line stands for the continuum model results 
with K « 0.39 (« = 0). The shaded region indicates the range 
of Poisson ratios from = 0) = 0.16 to = 1/2. 


The strain-displacement relations presented above sig¬ 
nificantly modify the electronic properties of deformed 
2d materials, which are crucial for strain-engineering [31 
1 HZ]. Essentially, mechanical deformations have two 
effects on the electrons. Firstly, the change in electron- 
ion potential in the neighborhood of an atom generates 
an on-site potential often referred to as the deformation 
potential. Secondly, changes in the distances between 
neighboring atoms modify the overlaps of the correspond¬ 
ing orbital wave-functions and thus the electronic struc¬ 
ture. 

These effects are accounted for in the nearest-neighbor 
hopping Hamiltonian which provides an accurate de¬ 
scription of the low-energy 7r-bands in graphene |31|51IZ], 
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namely, 




^ = E 


3i;o(ao) + — E 

/In ^ 


^ik ‘ \ -j* 


^0 


-to(ao) E ( ^ “ 


<j,*> 


/3 j. 


Oo Oo 


■ ( 11 ) 


Here /3/ao = -io(ao)/io(oo), ff/ao = '^ 0 ( 00 ), where to(r) 
is the (distance dependent) hopping amplitude and vo(r) 
is the electron-ion potential. In the following we omit the 
overall on-site energy 3no(ao)- 

To obtain analytical insight, it is customary HEU to 
consider the low-energy limit of the tight-binding model 
above. By expressing Eq. (1^ in reciprocal space and 
taking its long wavelength limit, that is, by expanding 
H to linear order in momentum around the K and K'- 
points Ellll, one obtains an effective Dirac Hamiltonian. 
For K we write 


Hk = vpa-■ {p - A) + VB^2 , ( 12 ) 

where vp = 3aoto/2 ~ 10^ m/s, cr are Pauli matrices, p 
is the momentum, A is the vector potential, and z;d is 
the scalar deformation potential. The latter are given by 


_ /3(1-^) .^ 

jtlx — ^xy 

ao 

_ / 3(1 - k) _ 

-^y ~ 2 ( 2 q \^yy ^xx) ; 

■ud = ^ [1 + 2Kp(h)h] {Uxx + Uyy) . 


(13a) 

(13b) 

(13c) 


Hence, we find that the strain-displacement relations 
renormalize A by a factor 1 — k with respect to the stan¬ 
dard elasticity-induced gauge theory [3 H [MS]- For 
graphene, the material parameters [14j give 1 — At « 3/5. 
(Using the approximation k « 2j//( 1 -f v) we obtain 
1 — K « {1 — v)/{l + v) = 2/3.) Based on a VFM 
that drastically overestimates the Poisson ratio, previ¬ 
ous works [51 ITU] report a reduction factor of about 1/3. 
Note that in the derivation above, we ignore the strain- 
dependent renormalization of the Fermi velocity reported 
previously [ 21 ], since it has a negligible effect on our re¬ 
sults. 

The reduced vector potential has an important effect 
on strain engineering applications [281130] . Let us analyze, 
for instance, the first (sub)band-gap of AGNRs[?7]. In 
the long wavelength limit, the magnitude of this band 
gap is given by 


E„ 


V3\to\ min 

n=0,l,... 


A2 + 1 
V3/1 


{n-4>) 


D 


(1 - k)(1 -I- v)Ux 


(14) 


with (f> = 0 for metallic and </ = 1/3 for semiconducting 
AGNRs and where N 2 is the number of unit cells in the 
transversal direction and n is the band index. The mag¬ 
nitude of the band-gap as a function of strain is piecewise 
linear and shows a characteristic zig-zag shape. We pre¬ 
dict that the slope of the linear regions is decreased by 
a factor (1 — «;)(« 3/5 for graphene) due to the strain- 
displacement relations, while the positions of the max¬ 
ima and minima of the band gap are shifted by a factor 
1/(1 — k)(« 5/3 for graphene). 

We verify this prediction numerically. First, we ob¬ 
tain the relaxed lattice structure of a stretched AGNR 
as described previously. Subsequently, we diagonalize 
the Hamiltonian 0 using to = 2.8 eV, g = 4 eV, 
p = 3.37[3i, and the relaxed lattice structure as an 
input. From the eigenvalues we directly calculate the 
band-gap shown in Fig. For comparison we also 
show Fg for K = 0. For the considered range of applied 
strain we find a good agreement between the numerical 
results and Eq. (14). We observed that the strain re¬ 
quired to achieve a certain gap size increases for k > 0 . 


A similar renormalization is also found for phospho- 
rene. Based on the two orbital tight-binding model 
put forward in Ref. [3 and the strain-displacement rela¬ 
tions derived in the supplementary material we find an 
anisotropic renormalization of the strain-induced band- 
gap. It is reduced by factors 2.5 and 1.35 when straining 
in X and y direction, respectively. 


III. CONCLUSIONS 

Two-dimensional materials have gained considerable 
interest, in particular due to the enticing prospect of en¬ 
gineering their electronic properties using strain. Our 
results show that the relationship between strain and 
electronic properties is less trivial than often assumed. 
For Bravais lattices, the relation between strain and dis¬ 
placements is simple, namely the bond vectors transform 
as r.y —)■ (U 3 -I- u) • Tij. For non-Bravais lattices, the 
basis atoms introduce additional degrees of freedom, sig¬ 
nificantly modifying the relation between strain and dis¬ 
placement. With few exceptions [51 ITT)] . this fact is typi¬ 
cally neglected in the literature on electronic properties 
of deformed 2d materials [51 HI 151171 E5H5D] . We show that 
under an applied strain the bond vectors of 2 d materials 
generally transform as (Ua -|- u) • -|- Ay -|- A±z. 

The vectors Ay and A±z are determined by minimiz¬ 
ing the deformation energy. For graphene, the strain- 
displacement effect on the electronic properties is to 
renormalize the vector potential by a factor (1 —k) « 3/5 
while keeping its functional form. We provide a simple 
estimate of this correction based on the macroscopically 
observable Poisson ratio p, as (1 — k) « (I — p)/{ 1 + v). 
The (I — At) renormalization alters the dependence of the 
band gap in AGNRs on strain: it increases the strain 
required to reach a certain band gap. Having a well 
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applied strain Uxx applied strain Uxx 

FIG. 4. Strain dependence of the (first) gap size for AGNRs of different widths {Ni = 300 in all cases). Symbols represent 
tight-binding calculations without (blue) and with relaxation (red). The dashed (solid) lines show the effective Dirac equation, 
Eq. ( [T4| ), results for k = 0 (k « 0.39). 


controlled band gap is of key importance in a variety 
of applications [21 HI [52H55] . 

For 2d materials with a more complex crystal structure 
than graphene, such as the TMDCs and phosphorene, the 
strain-displacement relations we put forward are key to 
understand their fundamental material properties. For 
phosphorene, we show (see Supplementary Material) that 
its anisotropic mechanical properties |3l [36] can only be 
obtained by considering proper strain-displacement rela¬ 
tions. Furthermore, phosphorene is considered to be very 
interesting from a strain-engineering perspective as it dis¬ 
plays a strain dependent band-gap [33139]. Our strain- 


displacement relations lead to an anisotropic renormal¬ 
ization of this gap. This Letter provides a transparent 
approach to estimate the electromechanical coupling in 
2d materials, given that the mechanical and electrical 
properties are reliably parameterized. 
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Supplemental Materials: Strain-displacement relations and strain engineering in 2d 

materials 


1 


STRAIN-DISPLACEMENT RELATIONS FOR BLACK PHOSPHOROUS 


Black phosphorous is a two-dimensional material with an orthorhombic lattice-structure with lattice vectors ai and 
a 2 , and four basis atoms arranged in a puckered structure (see Fig. SI). The bond vectors have approximately equal 
length d « 2.24 A, and intra- and inter-pucker angles ©i = ©328 ~ 96.3° and ©2 = ©123 ~ 102° [ST]. Hence, the 
bond vectors read 


ri 2 = [-dcos(© 2 )/cos(©i/2),0,d-\/l - (cos(© 2 )/cos(©i/2))^] , (Sla) 

r 23 = [dcos(©i/2),dsin(©i/2),0] , (Sib) 

r 43 = [dcos(©i/2), —dsin(©i/2), 0] . (Sic) 

with r 43 = r 28 = rgi, r 23 = r 4 i = ryi, ri 2 = r 54 and rsg = [ri 2 ,x,ri 2 ^y, —ri 2 , 2 ]. A lattice deformation, cast in terms 
of the strain tensor u, modihes the bond vectors as 


ri2 = ri2 + u • ri2 -f A|| -h Aj_i , 

^23 = I’23 + U- r23 - A|| , 

^43 = r43 + U • r43 - A|| . 


(S2a) 

(S2b) 

(S2c) 


We assume that = r'gi, r '23 = r^i = = r '54 and r(jg = -r'i 2 ,z\- 

According to the VFM put forward in Ref. [S2| (see also Ref. IS1|1 the elastic energy of the unit cell is given by 


F^cell = Krd^ {5r23f + Krd^ {Sr28f + K'rd^ {Sr2if 

+ 2KQd^ (d©328)^ + 2KQd^ (d©32i)^ -f 2KQd^ (i5©82i)^ 

+ 2Krr'd^ [5r23dr28) + 2K'„,d^ (drggdrgi)-I-2A:A/d^ {Sr28Sr2i) 

+ 2Kred^ (dr23d©234) + 2Kred^ (dr43(5©234) + 2K'^Qd^ (dr23d©i23) + 2K’^Qd? (dr28d©i28) 

+ 2K'^Qd^ ((5ri2(5©i23) + 2Ar(,0d^ ((5r’i2(5©i28) , (S3) 



FIG. SI. Lattice structures of phosphorene. It is characterized by two lattice vectors and four basis atoms (1,2,3,4), defining 
ten nearest neighbor vectors. 
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V. 

Yy 

Gxy 

without correction 

58 

95 

51 

with correction 

17 

94 

20 


'^XX 

"^yy 

Vxy 

without correction 6993 

8930 

5957 

with correction 

3508 

8147 3707 

Kaneta et al. |S2| 

3720 8190 

3760 


TABLE SI. Young’s and shear moduli and sound velocities from Eqs. (S6| and (S7|. The former are given in units of Njm 
and the latter in m/s. 


where ivij = |rL — Vij\/d « (rL — Vij) ■ is the relative bond-length change and 5Qijk ~ ~(cos(0'^-^) — 


- ^^Md « (rb - • Y,j/(P 

cos(0yfe))/sin(0ijfe) withcos(0b^) « {r'^j-r'ji.){l—Srij/d—6rjk/d)/(P is the bond angle change due to the deformation 
up to first order in strain and A. There are nine force-field parameters, namely, K'^, Kq, Kq, K„'^ Kre, K'tQ 
K"q- Those determine the energy cost for bond stretching, angle bending, bond-bon d a nd bond-angle correlations. 
The expressions for the bond vectors, Eqs. (SI) and (S2), render the energy Eceii, Eq. (S31, to be a quadratic function 
in the components of the strain tensor and the components of A. By minimizing Eceii with respect to Ay 2 ,, Ay j, and 
Ax and dividing by the area of the unit cell, Aceii, we obtain the deformation energy-density fceii = Eceii/^ceii- For 
the VFM parameters reported in Ref. we find 


Ay _a; = 0.61d Uxx + 0.37d Uyy , Ay = 1.02d Uxy , Ax = —0.36d Uxx — 0.20d Uyy . 


(S4) 


These expressions, together with Eqs. (S2|, constitute the strain-displacement relations for phosphorene. 


A. Elasticity constants 


Having the deformation energy-density fceii one can calculate the elastic constants Cij, which are given bv[S3| 


, _ d'^Scell 

duiduj 


(S5) 


where i,j = 1,2,6 and ui = Uxx, U 2 = Uxx and uq = Uxy From those one derives Young’s and shear moduli, Yxjy 
and Gxy^ and the Poisson ratios Vxyjyx-, 


Yx = 


IIO 22 — 


c, 


22 




/^2 

11O22 — ^12 


c 


11 


Gxy — Gee 


G 22 


'yx 


Cl2 

Wi ■ 


Moreover, the sound velocities are given by[S4j 


'^xx — 




'^xy — 



(S6) 


(S7) 


where pp is the 2d mass density of phosphorene. 

The results we calculate for all those elastic properties are summarized in Tab. |S1[ We also provide results for the 
case where the correction A is not taken into account. Comparing the sound velocities to the values given in Ref. 
M one readily concludes that is necessary to account for A to get comparable results. One may also compare the 
Young’s moduli to recent DFT results from various groups. There the values of the reported Young’s moduli are for 
instance: 26 N/m and 88 N/m [S5| . 24 N/m and 92 N/m [S3], 29 N/m and 102 N/m [S6| . Again the values in Tab. 
[STjfor the corrected case match best. 

Qualitatively, the VFM reported in Ref. |S2| gives very good results. However, it should also be mentioned that it 
has some deficits. For example, the Poisson ratios are obtained as i^yx = 0.27 and I'xy = 0.05, whereas other studies 
report 0.81/0.24[S5] or 0.62/0.17[S3]. However, it should be kept in mind that the parameters of the VFM were found 
by fitting to optical F-phonons measured in an experiment. To get a better quantitative description more data is 
required and the Poisson ratios need to be included in the fitting procedure. 
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B. Electronic band-gap 

As argued in the main text, the strain-displacement relations are very important for correctly estimating the 
influence of strain on the electronic structure. To demonstrate this for phosphorene, we follow Ref. [S7] where a two- 
orbital tight-binding model was used to calculate the electronic band-gap of phosphorene. The band-gap energy is 
then given by 


^'gap — 2(ti + t2 + tg) , 


(S 8 ) 


where ti, t 2 and tg are the hopping parameters between atoms 2 and 3, 2 and 1 and 2 and 8 , respectively. For 
the undeformed lattice the values are ^ 1=^3 = —0.797 eV and = 2.393 eV. Assuming a distance dependence of 
the hopping parameters as t{r) oc 1 /r^, in accordance with the nature of the orbitals responsible for the relevant 
electronic bands, allows it to calculate the strain-induced modification of the band-gap AFlgap. To this end we use 
the strain-displacement relations, Eq. ([S2| and (p3, for the respective bond-lengths. We obtain 


without correction: AEgap « 1.899uxx + 3.538uyy , 

with correction: « 0 . 755 ^ 2 , 2 , -I- 2.612uyy . 


(S9) 

(SIO) 


For comparison, in Ref. IS7lthe band-gap modification AAgap « 1.8Q8uxx -I-3.507mj,j, was obtained by transforming all 
bond-vectors using the strain tensor only. Similar to the case of graphene, one obtains a non-negligible renormalization 
due to the non-Bravais nature of the lattice. 
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